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One of the severe problems encountered in high-power-solid-state laser systems is the thermal 
damage to laser rods and optical elements. One such type of damage is thought to arise from metallic 
or dielectric inclusions; i.e., impurities with physical and optical properties which differ substantially 
from those of the host material. Such inclusions may absorb an appreciable amount of the incident 
radiation and thereby may undergo thermal expansion. This produces major stresses within the host 
material. Estimating such thermal properties requires the consideration of solutions to the heat dif- 
fusion equation and to the thermal stress equations with appropriate boundary conditions. The optical 
path length change for a probing light ray passing near the inclusion, the radial and tangential stress 
components, and the changes of the refractive index for radially polarized and tangentially polarized 
light due to the thermal stress field are computed. The dependence of the maximum value of the tensile 
stress upon the size of the inclusion and upon the physical properties of the host is examined. The 
feasibility of using optical techniques to detect metallic and dielectric inclusions in laser materials 
before they cause damage also is studied. The computations suggest that the use of laser pulse widths 
of the order of microseconds or longer may be more promising for the detection of small incipient 
absorbing centers than the use of nanosecond pulse widths. 
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1. Introduction 

One of the severe problems encountered in high- 
power-solid-state laser systems is the internal damage 
of the laser rod and its optical elements. One type of 
such damage is the formation of interior cracks which 
run more or less normal to the direction of the laser 
beam [l]. 1 This internal cracking is thought to arise 
from solid metallic or dielectric inclusions; i.e., im- 
purities with physical and optical properties which 
differ substantially from those of the host material. 
There are a few documentary cases for which the laser 
glass around inclusions had melted and other cases 
for which metal was found near the center of failure 
cracks [2, 3]. The metal has been usually platinum or 
antimony. The present theoretical study estimates the 
thermal stresses which occur when a solid inclusion is 
heated by the laser beam and examines the feasibility 
of using optical techniques to detect such inclusions 
before they cause damage. Detecting an inclusion 
before it causes damage eases the burden of the 
chemical analysis near the incipient damage area. 

E. Bliss fl] has outlined the general mechanism 
for internal cracking. A solid inclusion with a high 
absorption coefficient is heated relative to the host 
material. The resulting thermal expansion of the inclu- 
sion may stress the surrounding host material suffi- 



1 Figures in brackets indicate the literature references at the end of this paper. 



ciently for damage to occur. In fact, the theoretical 
model presented in this paper predicts that the heat- 
ing may produce stresses larger than the theoretical 
breaking stress of the glass. 

Most experimental data on inclusion damage has 
been obtained from laser glasses which have hot 
platinum near the molten glass, which are processed 
in platinum or ceramic crucibles, and which contain 
neodymium as the lasingion. The metallic inclusions in 
laser glasses can originate from many sources and can 
include most metals. However, considering the ther- 
mal properties of most metals that could occur as 
metallic inclusions in laser glass, one finds that 
they are similar in their capability for causing failure. 
The noble metals, in particular platinum or platinum 
alloys, are the metals that have been identified with 
any reasonable frequency in damaged laser glasses 
[4]. The nonmetallic inclusions may be divided into 
crystalline types and noncrystalline types, occluded 
or exsolved gas bubbles, and glassy regions. Non- 
metallic inclusions can produce several weakening 
effects in the laser host. The most apparent of these 
are the lens effect, strains in flaws, electrostriction and 
absorption. Few definitive statements can be made 
about these effects and the importance of each is 
evaluated on a probability basis [5J. However, most 
researchers expect that those metallic and dielectric 
inclusions with high optical absorption at the lasing 
wavelength are most likely responsible for failure 
of laser glasses. The numerical results of this paper 
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will treat the cases of platinum, antimony, and AI2O3 
inclusions as representing the thermal, elastic, and 
optical properties of several possible inclusions in 
laser glass hosts. The theory will be valid, however, 
for other hosts which exhibit isotropic elastic prop- 
erties and for other absorbing centers. 

The several properties of the center-host system 
which determine the probability for internal cracking 
to occur may be divided into four groups. The first 
group consists of the bulk properties of both the in- 
clusion and host and includes the respective thermal 
conductivities, heat capacities, thermal expansion 
coefficients, elastic properties, absorption coefficients 
for the incident radiation, and equations of state for 
the liquid and gas phases. The second group contains 
geometric properties of the inclusion which also 
influence its ability to cause fracture, such as the size, 
shape, and orientation to the incident radiation. The 
distribution and nature of initial microcracks and of 
optical imperfections form the third group of prop- 
erties which determine the resistance of the host to 
internal cracking. Finally, the fourth group of prop- 
erties describes the absorbing center-host interface; 
namely, the absorptance, emissivity, and initial ther- 
mal contact between the absorbing center and the 
host. The elastic properties, the thermal properties, 
and the optical properties of Pt, Sb, ALO.j, and two 
representative neodymium doped laser glasses are 
cited in tables 1,2,3,4, and 5. 



Table 3 

Values of parameters pertaining to neglecting the coupling and 
inertia terms of eqs (1) and (2). The quantities v, a 2 , {8IT ), and t r 
are respectively the speed of elastic waves, diffusivity, coupling in 
units of the ambient temperature 7o, and the relaxation time for the 
definition of a local temperature to be meaningful. 





V 


a 2 


(Wo) 


t r 




(cm/s) x 10 5 


(cm 2 /s) 


(1/°C) x 10-' 


(s) x 10- 13 


Pt 


3.67 


0.238 


4.39 


2.0 


Sb 


3.61 


.133 


1.67 


2.3 


Al 2 3 (a) 


11.44 


.0875 


4.83 


-10.0 


Glass(B) 


5.43 


.005 


1.11 


-10.0 


Glass(U) 


5.06 


.003 


0.417 


-10.0 



Table 4 

Optical properties for r« > \ ~ 1.06 jjl. The quantities A(k, 7V), 
n, n' , and y~ x are respectively the absorptance, the real part of the 
refractive index, the imaginary part of the refractive index, and the 
inverse of the absorption coefficient. 





A(k,To) 


n 


n 


(cm) X 10 " 6 


Pt 

Sb 

Al 2 0»(a) 

Glass(B) 


0.30 
.35 
.60 


1.14 

3.4 

1.67 

1.52 

1.52 


3.25 

4.1 

-0.0 

-.0 

-.0 


2.4 

1.9 
(see text) 
- 4 X 10 8 


Glass(U) 




- 4 X 10 s 









Table 1 

Values of the elastic parameters. The quantities E, G, v, and x 
are respectively Young's modulus, shear modulus, Poisson ratio, and 
isothermal compressibility. 





E 


G 


V 


X 




(N/m 2 ) X 10 10 


(N/m 2 )Xl0 10 




(m 2 /N)XlO- 10 


Pt 


14.7 


5.28 


0.39 


0.0449 


Sb 


7.78 


3.25 


.20 


.234 


Al 2 0,(a) 


41.4 


16.2 


.275 


.0326 


Glass(B) 


6.74 


2.79 


.218 


.249 


Glass(U) 


6.74 


2.79 


.218 


.249 



Table 2 

Values of the thermal properties. The quantities p, C, K, and a 
are respectively the density, the specific heat at constant volume, 
the thermal conductivity, and the linear thermal expansion coefficient. 





P 


C 


K 


a 




(g/cm :i ) 


(J/g °C) 


(W/cm °C) 


(l/°C)xl0-« 


Pt(B) 


21.73 


0.17 


0.67 


8.6 


Pt(U) 


21.73 


.13 


.67 


9.0 


Sb 


6.62 


.21 


.18 


11.1 


Al 2 3 (a) 


4.0 


.84 


.29 


11.0 


Glass(B) 


2.6 


.63 


.0084 


9.0 


Glass(U) 


3.0 


1.3 


.013 


9.0 



Table 5 

Values for the photoelastic coefficients, j)\\ and/>i 2 ,the stress-optic 
coefficients, B\\ and B ± . and the change of index of refraction with 
respect to temperature, (dn/dT). 





P11 


Pvi 


(m 2 /N) 
X10-' 5 


B L 

(m 2 /N) 
X 10- 1 " 1 


(dn/dT) 

(1/°C) 
x 10-* 


Glass(B) 
Glass(U) 


0.134 
.134 


0.225 
.225 


5.33 
5.33 


21.8 
21.8 


4.0 
4.0 



The model formulated in this paper contains many 
physical assumptions which are necessary to render 
the problem solvable. The major assumptions are 
summarized here and are discussed in greater detail 
in the following sections. 

(a) The inclusion is a sphere of radius r () and is 
always in good thermal contact with the host. The 
number of inclusions per unit volume is assumed to be 
sufficiently small so that they do not interact with one 
another. The effects of shape and orientation to the 
incident radiation also are neglected in the model. 

(b) The host material is isotropic, continuous, and of 
infinite extent. It also is initially at an ambient tempera- 
ture To and free from all stresses and strains. Because 
the energy content of the incident radiation is finite, 
the latter statement requires the temperature to be 
To at infinity and all stresses and strains to vanish 
at infinity. The distribution and nature of microcracks 
and optical imperfections are not treated in the model. 
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(c) The linear-thermal-elastic equations are assumed 
to give a reasonable description of the processes which 
ultimately may lead to catastrophic damage. These are 
coupled equations relating the temperature and the 
displacement vector from which the stresses and 
strains are computed. They are valid only when a local 
temperature exists and when distances are larger 
than atomic dimensions (~ 10 s cm). A relaxation 
time t r for the definition of a local temperature T is 
approximately the reciprocal of a characteristic vi- 
bration frequency of the material. These relaxation 
times t r for Pt, Sb, AI2O3, and the laser glasses are 
about 10~ ,;5 s to 10~ 12 s. Hence, the equations are mean- 
ingful only when times t are much larger than 10~ 12 s. 

(d) It is assumed that the radiation of heat by the 
center-host interface and by the heated glass close to 
the absorbing center may be neglected in the thermal- 
elastic equations. A black body at 600 °C produces an 
energy flux of 0.735 (J/cm-s). The laser beams for the 
examples studied in this paper contain energy fluxes 
at least 10 4 times greater than 0.735 (J/cm-s). The 
temperature gradients which occur from the time of 
cessation of the laser pulse to time of maximum optical 
path length change are always sufficiently large so that 
the energy flux due to thermal conduction is greater 
than the energy flux due to radiation. The present 
calculations show that the energy flux due to thermal 
conduction becomes small for times greater than a 
second and that the radiation effect should be included 
in the thermal-elastic equations whenever t > 1 s. 
Because the temperature is close to the ambient tem- 
perature whenever the time is greater than a second, 
the long time behavior is not in the region of practical 
interest for detecting the incipient damage centers 
before they cause damage. Hence, it is assumed that 
all times are less than a second. 

(e) The linear-thermal-elastic equations contain a 
coupling term and an inertia term. The coupling term 
and the inertia term may be neglected whenever the 
three characteristic times which occur in the absorbing 
center-host system satisfy a set of inequalities. These 
times are the following. The pulse width r of the inci- 
dent radiation determines in part the rapidity of heat 
generation. The characteristic relaxation time for 
temperature equilibration (thermal diffusion) is £7, 
which is of order (r-jd 1 ), where r is the radial distance 
and a 2 is the thermal diffusivity. The characteristic 
mechanical time required for the production of stress 
waves is t u which is of order (r/v), where v is the speed 
of propagation of elastic waves. The relationships 
among the coupling term, the inertia term, the time 
history of the displacement vector, and the time history 
of the temperature are complex. The easiest way in 
which to study these relationships is to solve exactly 
the linear-thermal-elastic equations for a simple one 
dimensional problem. Boley and Weiner have solved 
such a problem [6|. They demonstrate that when 
r > tu and tr > £a/, then the coupling and inertia terms 
are small compared to the other terms in the equations. 
The data cited in table 3 for the inclusions and hosts 
discussed here reduce the above inequalities to in- 
equalities containing the pulse width r and radial 



distance r. The data show a maximum value for 
£.1/ ~ 2.8X 10~ H (s/cm) r and a minimum value for 
t T ~ 4.2 (s/cm 2 ) r 2 . Hence, these data lead to the in- 
equalities containing r and r; namely r > 2.8 X 10~ H 
(s/cm) r and r> 6.7 X 10~ 7 cm. It is assumed that these 
inequalities are satisfied and thereby that the coupling 
and inertia terms may be neglected. The results to be 
discussed in this paper are based upon a linear theory 
and upon a model whose elastic, thermal, and optical 
properties are constants. The values for these prop- 
erties, cited in tables 1 through 5, are valid for tempera- 
tures near 300 K and for small strains. This treatment 
is not expected to be correct near the onset of damage. 
The problem is certainly a nonlinear one near the 
region of catastrophic damage. Whenever any of the 
assumptions become invalid, then the results should be 
viewed as suggesting trends in the behavior of the 
system. The inclusions which occur in laser glasses 
are most likely not the spheres for which the model has 
been formulated. They could be irregularly shaped 
globules. Some Pt inclusions occur as hexagonal 
platelets. Even though the model does not take into 
account such geometrical aspects, one hopes that it 
does give a reasonable description of the actual system. 
One also hopes that the manner in which it suggests 
one should proceed to increase the damage threshold, 
though probably not quantitatively correct, is qualita- 
tively correct. 

In section 2, the linear-uncoupled-quasi-static thermal 
elastic equations are applied to absorbing inclusions 
in laser materials. The model for studying spherical 
inclusions is then formulated in section 3 within the 
context of the quasi-static thermal elastic theory. The 
heat diffusion equation with time dependent boundary 
conditions for the temperature as a function of time t 
and of radial distance ris solved by Laplace transforms. 
Expressions for the optical path length change for a 
probing beam of radiation passing near the inclusion, 
the radial and tangential components of the stress 
tensor, and the stress induced changes, of the refractive 
index for radially and tangentially polarized light are 
developed in section 4. The numerical results and con- 
clusions are presented in section 5. Section 5 also con- 
tains an estimate for the lens effect of the heated region 
surrounding the absorbing center. Finally, those optical 
properties which determine the absorptance and emis- 
sivity of the inclusion are discussed in the appendix. 



2. Thermal Elastic Theory 

Consider an absorbing sphere of radius r {) imbedded 
in an initially isotropic continuous host of infinite ex- 
tent. The thermal, elastic, and mechanical properties 
which enter linear thermal elastic studies are: the 
density p(g/cm ; *), the specific heat C at constant volume 
(J/g °C), the thermal conductivity K(W/cm °C), the 
linear coefficient of thermal expansion a(l/ °C), and 
the isothermal Lame elastic constants, A. and /x. The 
Young's modulus £(N/m 2 ), the shear modulus f/(N/m 2 ), 
the compressibility x( m2 /N), and the Poisson ratio v 
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are related to the Lame constants X and fx as follows, 
E= (3X + 2/x)/Lt/(X + jLt), G = £/{2(l + i/)}, 



X = £" 1 =3/(3X + 2 ) u), and 



:X/{2(X+*,)}, 



where the bulk modulus is A(N/m 2 ). The derived quan- 
tities of thermal diffusivity a 2 = (K/pC) (cm 2 Is) , the 
volume specific heat pC(J/cm 3 °C), and the speed of 
propagation of dilatational waves in an elastic medium 
v — '{(X -h 2p,)/p} 1/2 also enter the linear thermal elastic 
theory. The subscript c on any of these thermal prop- 
erties will refer to the absorbing center and the sub- 
script h will refer to the host. 

The basic theory which describes the behavior of 
an absorbing sphere imbedded in a continuous medium 
and subjected to a laser beam combines the theories 
of elasticity and of heat conduction under transient 
conditions and is a boundary value problem of con- 
siderable mathematical difficulty. Fortunately, 
simplifying assumptions without significant error 
are possible. 

The problem is that of one absorbing spherical in- 
clusion imbedded in an infinite medium initially at a 
uniform temperature 7o and subjected to a prescribed 
rate of internal energy absorption per unit volume 
Q v (r, t)(W/cm 3 ). The quantity Q v is related to the 
intensity of the local radiation / by an absorption 
equation; that is, 

Q v (r,t) = yl(v,t), 

where the linear absorption coefficient is y in units of 
inverse centimeters. This y should include also the 
effects of nonlinear absorption processes which be- 
come important at high intensities. The local intensity 
for radiation propagating in the plus z direction at z 
is related to the local intensity at Zq by the relation 

I(z) = exp {— y(z — Zo)}I(zq), 

where z > Zq. 

The linear-thermal-elastic equations relate the 
temperature T(r, t) and the displacement vector 
u(r, t) as functions of position r and time t. These 
coupled equations are [6] 



KV 2 T - pC ^ - (3X + 2p,)aT, j t (V • u) 



and 



d 2 u 



+ Qv = 0, (1) 



j*V 2 u + (X + /x)V(V-u)- P — 

- (3X + 2p,)aVr=0, (2) 

where To is the reference temperature for which the 
material is free of stress. The reference temperature is 
chosen to be 293 K. The strain and stress components 



are obtained from the displacement vector u. Equa- 
tions (1) and (2) are meaningful only when a local 
temperature exists. 

A close relationship exists between neglecting the 
coupling term (3X + 2fJi)aTo[d(V • u)/dt] and ne- 
glecting the inertia term p(d 2 u/dt 2 ) in eqs (1) and (2). 
This relationship is based upon the physical idea that 
the rate of thermal expansion may not be sufficient to 
produce elastic waves. Let us consider the coupled 
heat equation. We rewrite eq (1) in the form 



KV 2 T - pC 



c]T 
dt 



1 + 8 



(H2/x){a(V>u)W 
(3k + 2fji){a(dTldt)} 



+ Qc = 0, 



d') 



where the dimensionless parameter, 



8={(3X + 2n) 2 a 2 T () lp 2 Cv 2 }, 

determines the amount of coupling. Equation (1 ') shows 
clearly that neglecting the coupling term requires that 
8 must satisfy the inequality 8 < 1 and that the strain 
rates, d(V • u)/d£, must be less than or of the same 
order of magnitude as temperature rates a(dT/dt). The 
latter condition implies that the time history of the 
displacements must follow closely that of the tempera- 
ture. In fact, the authors of ref. [6| have shown that if 
the application of heat is not too rapid, not only may 
the coupling term be omitted but the inertia term may 
be omitted under the same circumstances. 

The rapidity of heat generation at the absorbing 
inclusion-host interface is characterized by the laser 
pulse width r (s). In addition to the time r, two other 
characteristic times exist in the system. Namely, the 
characteristic thermal time, tj— (r 2 /a 2 ), is the relaxa- 
tion time for temperature equilibrium (thermal diffu- 
sion) and the characteristic mechanical time, t\i= (r/v), 
is the time required for the production of stress waves, 
where the quantity r= |r|. It has been demonstrated 
[6] that when 



r > tM 



and 



t T >t M , 



(3) 



then coupling and inertia terms are small. That is, when 
inequalities (3) are satisfied, then the linear-uncoupled- 
quasistatic thermal elastic theory obtains; namely, 



K57*T-pC^ + Q v = 0, 

Ot 



and 



p,V 2 u-h (X + /i)V(V • u) - (3X-h2/>i)aVr=0, 



(4) 



(5) 



where V 2 is the Laplacian operator, and V is the 
gradient operator. 

Hence neglecting the coupling term one finds that 
the thermal stress problem becomes two distinct prob- 
lems which are solved consecutively. The first problem 
is the solution to the boundary value problem whose 
field equation is given by eq (4). When we have the 
temperature distribution T(r, t) , the second problem 
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is the determination of the resulting stress distribution across the center-host interface is 
from eq (5), and the stress-strain relations, 



<jjj = 8ij\€kk + 2/Lt€ij — 8jj (3X + 2fx)aT, 



(6) 



where 6\-j = when i 4= j and 8n — 1 when i=j. Equation 
(6) is subject to the equilibrium equations, 

V<r,i + f=0 



€ij= (i)[(dui/dxj) + (dujldxi)], 

are the strain displacement relations, crtj is the stress 
tensor, and f is the applied force per unit volume. The 
subscripts i 9 j 9 and k represent the cartesian compo- 
nents x, y, and z. 



3. Spherical Inclusions 

The heat diffusion eq (4) describes the time and 
space dependence of the temperature 7"(r, /) for 
large distances r and large times t; 



a 2 V 2 T=(dT/dt). 



(7) 



Equation (7) does not contain the volume heat genera- 
tion term Q v (r, t) because the host is assumed to 
be very transparent to the laser energy and to absorb 
only a negligible amount of energy from the laser beam. 
In addition, because the absorption coefficient for the 
absorbing center y r (Pt or Sb) is very large compared 
to the absorption coefficient for the host y/, , y c ~ 10 ,4 y/,, 
the heat absorbed by the absorbing center from the 
laser beam is assumed to occur entirely at the center- 
host interface. This is not the case for Al 2 ;i because 
y r (Al 2 3 ) ~ Jh (glass). If the Al 2 0. r glass interface 
were highly absorbing due to some mechanism (Max- 
well Demon), then the assumption that heating occurs 
only at the interface would be valid. The nature and 
existence of such a mechanism for AI2O3 is specula- 
tive at present. Hence the author uses in this paper 
the name Al 2 0.*(a) to denote an inclusion whose 
properties correspond to those for polycrystaline 
ALO3 with the one exception that the absorption 
coefficient y r (Al 2 ;i (a)) is very much larger than 
that for y r (Al 2 O:0 ~ y/, (glass). 

The diffusion equation requires a statement of the 
boundary conditions before a solution is uniquely de- 
fined. The temperature has the form 

I To, t^O 

e(r,o= (8) 

\To + T(r 9 t), t>0, 

for all values of r. The function T(r, t) is the tempera- 
ture relative to the ambient temperature To. The 
boundary conditions are time dependent for any 
realistic treatment of the problem. The temperature 
must be finite everywhere; 7\r=o°, t) = finite and 
jT(r = 0, t) = finite. The continuity of the temperature 



T c (rof, t) =Th(ror, t) 



(9) 



where f is the radial unit vector. Condition (9) implies 
that the host is in perfect thermal contact with the 
inclusion for all times t 3? 0. The conservation of heat 
flow at the interface gives 



KcVTc(r„r, t) ■ t + H{t\ r • z] 



= — K f ,VT, l (n ) r, t) • r, 



(10) 



where z is the unit vector in the +z direction. 

The energy flux H(t; rz)(W/cm 2 ) is the energy 
absorbed by the sphere from the laser beam. Boundary 
condition (10) describes an uneven heating of a sphere 
by a plane wave propagating in the positive z direction. 
It contains angular properties which are unnecessary 
refinements for our purposes, particularly when the 
thermal conductivity of the center is much greater 
than that of the host. Therefore, it is reasonable to 
introduce one more assumption which reduces bound- 
ary condition (10) to a spherically symmetric one. Let 
El be the energy density (J/cm 2 ) of the laser beam and 
let t(s) be the pulse width of the laser beam. A square 
wave pulse is assumed for our problem. The spherical 
inclusion intercepts an area irrl from the laser beam 
and absorbs the amount of energy per unit time 
Q T =[A(K T)(EijT)]mi in watts, where A(K T) is 
the absorptance of the center. If the sphere were to 
absorb energy at the rate Q T over its entire area, then 
the equivalent energy flux //(W/cnr) would be 
H = ((?t/477To). The equivalent energy flux for the 
spherically symmetric case is given by the relation, 



H(t) = 







for 



«?r/47rr§) 



;<0 
for 



and 



t >r 



Hence, replacing the angular dependent energy flux 
H{t\ f'z)> by the spherically symmetric energy flux 
H(t) reduces the problem to a spherically symmetric 
one. Boundary condition (10) now reads 



■KcVT c (r n r,t) • r + H(t)=-K h S7T h (r K t) 



(11) 



Boundary condition (11) corresponds to replacing the 
incident plane wave, which has an energy flux (EJt), 
with an incoming spherical wave, which has an energy 
flux (7rr?£ L /47rr 2 r). 

The author does not include the radiation of heat 
by the interface and by the heated glass close to the 
inclusion in this model. The temperature gradients 
which occur from the time of cessation of the laser 
pulse to time of maximum optical path length change 
are always sufficiently large so that \Ki,VT},\ is greater 
than the energy flux due to radiation. Of course, the 
energy flux \KhVTh\ becomes small for times greater 
than a second and the radiation effect should be in- 
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eluded in both the diffusion equation and the boundary 
conditions whenever t > Is. Because the temperature 
is close to the ambient temperature whenever the 
time is greater than a second, the long time behavior 
is not in the region of practical interest for detecting 
the incipient damage centers before they cause 
damage. 

The solutions for the temperature T(r, t) are ob- 
tained by taking the Laplace transform of the diffusion 
eq (7), the initial condition (8), the continuity eq (9), 
and the time dependent boundary condition (11). The 
transformed equations become linear equations for 
the Laplace transform [/(r, s) of the temperature 
T(r, t), where s is the Laplace transform variable 
for the time t. In a separate paper [7], it has been 
shown that 

U h (r, s) = B fl exp (— z h )lz h , (12) 

where zn = (rs 1/2/a /,), 

B ft = h(s)(z h0 /s 1 l 2 ) exp (z h0 ) X [(Khloh)(l+Zkt) 

— (Kcla c ) (zeb 1 — ctnh z c o)]~\ 



h(s) = [A(K, r )£/74T][l-exp (-st)]*" 1 , 



Zh 



o = (r () s ll2 la fl ) and z c o = (r s ll2 la c ). 



(13) 



Evaluating the Bromwich contour integral to obtain 
the temperature F(r, t) from a knowledge of the 
Laplace transform £/(r, s) , eqs (12) and (13), would be a 
most ambitious task. 

Two different arguments exist by which eq (13) re- 
duces to a more tractable form. The first argument 
treats eq (13) in the limit that z c o< 1, then the ctnh 
(zco) in eq (13) is approximated by 



ctnh (z c o) -Zco 1 + Uco/3) +Q(z$o) + . . ., 
and eq (12) thereby is approximated by, 



Ui l00 (r, s)=Lim U,,(r, s) = J3/ t (o°){exp (-z„)/z fl }, 



(14) 



where 



B, ,(<*) = 



h(s)zho exp (zho) 



M{s l l 2 + b + )(s l l 2 + b-) ' 
M=(p c C c ro/3),6 + = 6[lT(l-/?) 1/2 ], (15) 

b= (Kkl2Ma h ) 9 and R = (4p c C c /3p/,C„). 

The second argument yields the same relations as eqs 
(14) and (15) but involves simplifying the boundary 
condition (10). Observe that because a 2 > ISaf and 
K c > 20Kh the temperature gradients in an absorbing 
center such as platinum, antimony, and Al-iO.^a) would 
be much less than those for the glass host. In the limit 



that K c > Ki,, the boundary condition (10) becomes 
Lim [M(dT c (ro, t)ldt)]=Lim K h (dT„(r+t)ldr) 

+ H(t), (16) 

where r± = ro±e. Hence, approximation (14) obtains 
either when \ros ll2 la c \ ^ 1 or when K c > K/,. The latter 
inequality imposes no restrictions on r () and 5, but the 
former inequality does impose such restrictions. The 
volume specific heat ratio R= (4p c C c /3p/,C/,) deter- 
mines three regions of behavior. Namely; when R < 1, 
the roots b ± are real and unequal; when R = 1 , the roots 
are real and equal; and when R > 1 the roots are com- 
plex and conjugate to one another. The evaluations of 
the Bromwich integral, 



1 f c + i(I 
T/,(r, t)=Lim- — : exp (st)Ut,(r, s)ds, 

rf-*« Irn Jc-id 



(17) 



for the above three regions appear in ref. [7], and give 
the temperature T(r, t). The quantity c is chosen to be 
sufficiently large in order that the integral, 



/: 



exp (— ct)\Ti,(r, t)\dt, 



(18) 



exists. Including the mathematical details for the 
evaluation of Ti,(r, t) in the present paper would only 
detract from the presentation of the model, of the 
physics which it contains, and of its numerical results. 
However, the temperature gradient in the host is nega- 
tive for all values of position and time and for all finite 
values of the volume-specific-heat ratio R. 

4. Theoretical Expressions 

Having thus obtained the temperature T(r, t), one 
may compute all the other quantities mentioned in sec- 
tion 1. The optical path length change AL becomes 



AL(n): 



■/. 



(dn h ldT h )T fl (r,t)dx, 



(19) 



where the integration path is along a ray that has a 
distance of closest approach n to the absorbing sphere 
and the variation of the refractive index of the host with 
temperature is (dnnldTh). The quantity (drthldTh) is 
properly the rate of change of the host refractive 
index with respect to temperature evaluated at zero 
strain. When (drthldTh) is independent of r and t, 
eq (19) has the form, 



AL(n; r {) ) = (dn h ldTh) f + ° T h ((r\ + x 2 yi 2 , t) dx, 

J-D 



(20) 



where T h ((f\ + D 2 ) 1 ' 2 , t) < T h (r , t) and the integra- 
tion path is a straight line of length 2D with a distance 
of closest approach to the sphere n ^ r . The stress 
birefringence has been neglected in eq (20). The 
quantities which enter the stress birefringence are 
computed at the end of this section. 
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The stress components satisfy second order differ- This temperature profile leads to the following expres- 
ential equations which follow by inserting the tempera- sions for u c , (Trr i a an d a"ee t c = 0"<w,c when r < r : 
ture from eq (17) into eq (5) and by using the relation 



(6). The displacement vector u= ur is 

B(r) _«o+£if' +c ^ 

(1 — v) r 2 Ja r 2 

where a is any fixed radius and where p is the radial 
variable of integration. Equation (21 ) is correct only for 
those cases in which [{a(\ + v)l( 1 .- v)} (dT/dp)] + 
for a ^ p ^ r. Because the total energy (Q t t) deposited 
by the laser beam into the center-host system is finite, 
conservation of energy gives 



and 



f/ r (r, t) = C u r 



CTrr,c(r, t) = cr^, r (r, t) = ov*>, r (r, t) 

E c 



(\-2v c ) 



{Cir — OLcThinu t)}, 



(28) 



(29) 



where the condition a r (°°, t) = has been used to 
determine C>> r — 0. The corresponding expressions for 
the region r ^ r are: 



(? T r=(47r/3)peC c rgn(ro,0+47rp^ f" T H (r n t)i*dr. U »^ = ^ -tjr 3 J ro r/l(p ' *)P* d P + ir> (30 > 

Whenever [{a(l + v)j{\ p)}(dT/dp)]=0 for ( 7 ,., %/ ,( r , =- T^g^ , P T h (p, t)p 2 dp 

a ^ p ^ r then the displacement vector u= ur becomes v* Vh)r J /() 



tf(r)=C,r+(Ct/r*). 



(23) 



When [{a(l + f)/(] -- e) }(d773/-)] =1= 0, then the 
radial stress component has the form 



. 2aE f £C 



2^) 



2EC, 
(]+f)r ! 



(24) 



and the tangential stress components tree and </>«« are 
equal and are given by the equation, 

'«('■) = (Ti^iJ. T P- d P + TV^T) 



+ 



fC, 



a£T 



(l+^r 3 (l-v) 



(25) 



When [{a(l + v)l(\ ■- v)}{dT/dr)]=0, then the 
following results, 



E 2EC> 



(l-2») 
and cree (r) =ov*>(r), 



(l + p)H 



(26) 



o-»»(r) 



(l-2v) 



(r;,-«7') + 



ECt 

(l + v)r ; 



(27) 



The Laplace transform given by eqs (14) and (15) 
corresponds to the following temperature profile: 



T(r, t) = 



7 1 /, Oo, t) for 



T H (r, t) 



for 



=£ r ■■ 
r > r„. 



r<> 



2EhC->i, 
(1 + ^)r 3 ' 



(31) 



(1 - ^)HJ r() 



+ 



EhCou ai>E h T h {r, t) 



(32) 



where the condition a/, (°°, =0 has been used to 
determine C\/, = 0. Observe that eqs (31) and (32) 
predict that 07,-, /,(<», 0=0 and (Tee^i 00 , t) = 0. 

The coefficients C\ c and C21, are evaluated by in- 
troducing an adhesive boundary at r = rof. An ad- 
hesive boundary maintains the thermal contact between 
the absorbing center and the host and obtains when 
the following conditions are valid: 

(a) The hydrostatic pressure, 

P ( = — 3 (o- r/ ., r + o-00. ,. + ovv\ r ) = — ov,-, r , in the 
center at r=r {) r equals the negative of the 
radial stress component in the host at r=r<>r; 
that is o> r , c (ro, t) =cr rr , /i (r , *)• 

(b) The tangential stress components are discon- 
tinuous across the interface. 

(c) The radial displacement vector ur is con- 
tinuous across the interface; that is, //<(r , t) 
= Uh(r , t). 

Conditions a and c then yield two simultaneous equa- 
tions for C\ r and Con. The solutions to this set of 
equations are: 



C'2h 



a v Tu(r {) , t) 



l + {2£„xc/3(l + *>„)}' 



(33) 



and 



C\e= (C 2 i,lr'i). 
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A physically intuitive thought experiment which 
leads to the same expression for C211 exists. The sphere 
of radius r is removed from the host and is heated to 
Th(ro, t) at zero pressure, Pi = 0. After heating, it has 
an initial volume, 

Vi=(47rl3)r* [l + 3a c T h (r 09 t)]. 

The heated and expanded sphere is forced back into 
the host at a final pressure P/=— ov r , ft(ro, t) . It then 
has a final volume 



where 



^=(477/3)7-3(1+36), 

e= (u/ ? (ro)/ro) = (C 2 ft/rg). 
Combining the above relations gives 

o-rr,h(ro,t)=-{2E h el(l + v h )}. 

Finally, the isothermal compressibility \c for the center 
is introduced in order to compute e, 

X c=-aiv)(dvidP) T . 

Assuming the isothermal compressibility to be a slowly 
varying function of pressure for Pf^P^ P\ yields the 
relations, 

X c{Pf-Pi)=-\K{V f IVi) 

= -\n [(l + 3€)/(l + 3a c n(r ,0)]- 

When 3e < 1 and Sa c Tn (r , t) < 1, the expansions of 
the logarithm give the strain e; namely, 

a c Tn(r , t) 



l + {2E*Xc/3(l + v*)} 



Substituting the coefficient C2h = r'ie and the tem- 
perature T/,(r, t) into eqs (31) and (32) gives the radial 
and tangential stress components. The radial rates of 
expansion during heating and of contraction during 
cooling for the center also are computed, 



_de_ 
r °dt = 



roa c {dTh(ro, t)/dt) 
'l + {2E hX clHl + Ph)} 



The latter rates are used to check that the linear-quasi- 
static-thermal elastic theory is predicting self consist- 
ent results. If the rates of expansion or contraction of 
the center were to approach or exceed the speed of 
sound in the host, then the above problem becomes a 
nonlinear one involving two coupled differential equa- 
tions. One equation treats heat diffusion and the other 
treats energy transfer in the form of sound waves. 

One goal of this study is to examine under what 
conditions the absorbing center may initiate damage 
in the host. When the tensile stresses produced by the 
heating exceed the strength of the host, then damage 
should occur. The theoretical breaking stress is esti- 
mated theoretically [8] to be about one tenth of the 



Young's modulus, o^/, ~ {Eh 110). This value does not 
take into account the presence of microflaws and there- 
fore is a representative upper limit for the strength of 
the host. The sign convention of the foregoing formal- 
ism is that compressive stresses are negative and 
tensile stresses are positive. The radial stresses are 
according to eq (31) always compressive (negative) and 
hence do not initiate damage. The tangential stresses 
are according to eq (32) either compressive or tensile 
and depend upon the time, the position, and the prop- 
erties of both the center and host in a most complicated 
manner. Whenever the tangential stresses are tensile 
(positive) and are sufficiently large, 

(Teo,h(r, t) kath ~ (En 110), 

they may initiate damage. 

The tangential stresses at the interface may be 
written in the form, 



0-00, /,(r , t) =EhTh(r , t)da en , 
where the effective expansion coefficient is 

OLc Oih 



5a P 



1 + *a+(2EaXc/3) (l-vh) 



(34) 



(35) 



It is interesting to observe that even when a c = a/,, 
o"06>, /<( r o, t) =4=0. The stress cree,h is usually propor- 
tional to (a c — OLh) for those problems in which the 
temperature is spatially uniform everywhere. But in 
the present problem, when r ^ r , the temperature is 
not spatially uniform (that is, when r ^ r , V7 7 /, =N0) 
and the stress dee, n is not necessarily proportional to 
(a c — a/,). The model by which one reasons intuitively 
that the tangential stress is proportional to (a c — ah) 
is based upon eqs (23), (26), and (27) for both r < r 
and r > r<>. This is not correct when (dT/dr) differs from 
zero for r ^ r<> as it does in the present model. 

Whenever (dT/dr) = for ^ r ^ 00, a different set of 
boundary conditions must be used to determine the 
constants in eqs (23), (26), and (27). To obtain the 
correct boundary conditions for the infinite host, one 
proceeds as follows. Consider a sphere of radius r 
inside a spherical shell of inner radius r () and outer 
radius r 3 . The system is heated adiabatically from 
To to To + T so that (dT/dr) = for ^ r ^ r 3 . The cor- 
rect equations are eqs (23), (26), and (27) and the cor- 
rect boundary conditions are u c (0) = 0, u c (ro) = u/,(r ), 
cr,r, c(ro) =or rr , /j(r ), and 07,-, /,(r 3 )=0. These equa- 
tions and boundary conditions then predict that in the 
limit as r 3 approaches infinity, the tangential stress 
approaches the value, 



cree,h(r; r- 3 = oc) = 



E h (a c — a h )Trl 
{\ + v h +(2E hX c/?>)}rV 



(36) 



Equation (36) is not correct whenever (dT/dr) differs 
from zero. Observe that the boundary conditions lead- 
ing to eq (36) replace Uh( co )=0 with o- rr ,/j( 00 ) = 
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and also yield M 00 ) = °° an d cree.hi 00 ) = 0. Compare 
this with the results from eqs (30), (31), and (32); 
namely, u/,(o°)=0, cr r r,/j( x ) = 0, and or^,r(°°) = 0, 
because T/, (r, approaches zero faster than (II r 2 ) 
approaches zero as r approaches infinity. 

The maximum tensile stress occurs, when (dcree.t (r, 



t)/dr) evaluated at 



is zero and when (d 2 cree,h{r, 



t)ldr 2 ) evaluated at r = r max is negative. These 
conditions do not necessarily obtain at r = r max = r 
and at t = T. The qualitative behavior of the tangential 
stress depends in part upon the sign of the effective 
expansion coefficient da eff . When 8a eff is positive and 
sufficiently large, 8a e ff > 10~ 6 (1/°C), then the maxi- 
mum tensile stress does occur at r = ro and t = r. 
For fixed times t, the tensile stress decreases mono- 
tonically to zero as r increases from ro. When 8a e{f is 
negative, the maximum tensile stress does not occur 
at r=r and it also does not occur necessarily at t = r. 
Instead, the tangential stress is always compressive 
at r = ro, when 8a ef , < 0. Keeping time fixed one finds 
that as r increases from r the tangential stress in- 
creases algebraically from negative (compressive) 
values to positive (tensile) values, attains a maximum 
tensile stress at r= r max > r () , and then begins to 
decrease algebraically. When 8a eff is small, |8a eff | 
< 10~ 6 (1/°C), then the stress may change from tensile 
to compressive or from compressive to tensile more 
than once as r increases from ro. This behavior for 
small 8a e{{ is very sensitive to the sign and value of 
(dcroojiir, t)/dr) evaluated at r = r . Because evaluat- 
ing T f ,(r, t) and its derivatives to study analytically 
the maxima of croejiir, t) would be a monumental 
task, the numerical results from a computer have 
been used to obtain the above qualitative features. 
The effective expansion coefficient Sa eff is about 
-5.5 X 10- 6 (1/°C) for Pt in Glass (B) and it is about 
-6.7X 10-«(1/°C) for Sb in Glass (B). 

Finally, the expressions for the strain induced change 
of the refractive index for radially polarized light An r 
and for tangentially polarized light Arte = An* (stress 
birefringence) are developed. The theory of ref. [9] 
which examines the stress birefringence of an isotropic- 
continuous transparent medium is used. The strain 
induced change of the refractive index for radially 
polarized light is 

A/i r (r, t) = B\\o- r r,f(r, t) -\-B + (ae8,h(r, t) 

-fov*>, ft'(r, t))+a ft T h (r, (pu + 2pi 2 ), 

and the strain induced change of the refractive index 
for tangentially polarized light is 

An 6 (r, t)=An<p(r, t) =B±<j n , /,(r, t) 

+#|| o-H8,h(r, t) +B±<r<p<p,n(r, t) 

+ a h T h (r, r)(pu + 2pi 2 ), 

where the stress optic coefficients are 

B\\=Eh 1 (pu—2v h pi2), 



B ± =Ej; 1 [(l — vn)p a — vnp ii\, 

and where the photoelastic coefficients are pw and 
P12; i.e., 

An r = n r — n(e = 0) =pne rr -hpi 2(6^ + 6^), 
and 

Ane = ne — n(e = 0) =p ]2 (e r ,- + €^) +Pn€$e. 

The stress birefringence at any point r with respect to 
the local unit vectors for spherical coordinates is 
(An,- - Arte). 

Computer programs to evaluate the quantities 
AL, (Trr, O"00? Ara,, and An^ as functions of r and t for 
platinum, antimony, and ALO.^a) spherical inclusions 
are developed. Two different host glasses with prop- 
erties such that both R < 1 and R > 1 are considered. 
The strain rate ro(de/dt) and the energy flux 
S fi (Tf,(r, t)—T$) radiated by the spherical surface 
and by the host near the inclusion also are computed. 
The Stefan-Boltzman constant is denoted by S#. 
For the values of radial distance and time which are 
important to the inclusion damage problem, the effects 
of both quantities are negligible and hence the model 
is numerically consistent with its assumptions. 

5. Results and Conclusions 

The numerical results predicted by the model devel- 
oped in section 3 are reported in this section. Among 
the many input parameters the absorptance A(k, T) 
is perhaps most sensitive to the initial thermal contact 
and surface conditions of the host and inclusion. To 
avoid questions about the value of A(k, T) the 
numerical results for the reference (normalized) case 
when A(\, T) = l are given. This presents no addi- 
tional problem because the temperature Ti, (r, t), 
the optical path length change AL, the stress com- 
ponents (Trr and (Tee—trw and the stress birefringence 
Arc, and An^, are all directly proportional to the 
absorptance. 

The model of section 3 has been used to study two 
different questions. First, how the maximum tangential- 
tensile stress varies as a function of the radius of the 
spherical inclusion for a fixed energy density and pulse 
width of the laser beam. Second, if the maximum 
temperature of the inclusion is limited to a fixed value, 
Ti,(r=r{), t = r) = constant, how the parameters of the 
incident laser beam should be varied to increase the 
probability of detecting by optical techniques a small 
incipient absorbing center, before it causes damage. 
Optical techniques to detect small inclusions become 
more promising the greater AL, An,, and Ane are. 

Table 6 for inclusions in Glass(B) and table 7 for 
inclusions in Glass(U) contain the variation of the 
maximum tensile stress cr^(max-tensile) as a function 
of r for a laser beam having an energy density 20(J/cm 2 ) 
and a pulse width 30 ns. When 10~ 4 cm > r > 5 X 10" 5 
cm the maximum tensile stress exceeds the theoretical 
strength of the Glass(B). Also when 10~ 4 cm > ro > 5 
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X 10~ 5 cm the maximum tensile stress exceeds the 
theoretical strength of the Glass(U). These results 
demonstrate that submicron-sized inclusions have the 
greatest probability to cause damage in laser glass 
hosts. The results given for r = 10 -7 cm may have no 
physical meaning because the uncoupled thermoelastic 
theory becomes suspect for r < 6.7 X 10 -7 cm. When- 
ever the temperature Ti,(r {) , t) exceeds the strain point 
of the glass, the results reported in tables 6 and 7 are 
at best order of magnitude estimates for the tensile 
stress. Laser glasses have strain points around 600 °C. 
These two tables show that very large and very small 
inclusions are not likely to produce damage. Only 
inclusions having an intermediate size of the order of 
tenths of microns have a high probability to produce 
damage. They also show that the surface temperature 
Ti,(r {) , t) and the maximum tensile stress become 
independent of the thermal and elastic properties of 
the inclusion when r < 10~ fi cm. This observation is 
proven in ref. [7] by analytic methods. 

The high temperatures which are predicted by this 
model for 0.1)H ^ r ^ 1/jl are certainly not to be taken 
seriously. The model does not treat the heats of fusion 
H F , the heats of vaporization flV, the changes in volume 
upon melting, and the equation of state for the liquid- 
gas transition at extremely high temperatures. These 
equations of state for platinum, antimony, and Al 2 0,j 
are not known. Also some researchers argue that these 
effects which the present model neglects tend to cancel 
one another in the case of platinum for temperatures 
up to its boiling point and that thereby, even though 
the model predicts temperatures which are too large, 
it does give representative values for the stresses. 
Neglecting the possible chemical and physical changes 
which may occur at the center-host interface, they 
argue as follows. The heat of fusion H F (c) gives an 
effective temperature of fusion T F (c) = H F (c)/C c by 
which the temperature Tt, (r , r), given by the model, 
is too large whenever Tu{c) < Tf,(r , r) < T B (c), 
where Tm(c) and T B (c) are respectively the melting 
and boiling points of the center. The values for platinum 
are 7>(Pt) ~ 800 °C, 7\,(Pt) ~ 1770 °C and T B (Pt) 
~ 3830 °C. The excess temperature T f (Pt) predicted 
by the present model produces an excess volume 
change (&V/V) ex - Sa c T F (Ft) ~ + 0.02. But the 
volume change for the platinum upon melting is also 
about + 0.02 [9]. Hence, the model should give correct 
order of magnitude estimates for platinum up to its 
boiling point. Such arguments fail for antimony 
because its (AV/V) upon melting is about— 0.014. 

The following mesh in the r and t space is used to 
obtain the values given in tables 6 to 9. The expressions 
are evaluated for five values of the time (3 ns, 30 ns, 
0.3/xs, 3/xs, and 30/xs) and for twenty values of r. The 
limits on r depend upon r . When r ^ O.ljU,, the limits 
are r ^ r ^ 5r ; when r = O.Olpt they are r ^ r 
^ 60r ; and when r = 0.001/1, they are r ^ r ^ 600r . 
The twenty values of r are equally spaced in a loga- 
rithmic manner between the lower and upper limits. 

The quantities AL(r, r , t), <J rr {r, t), cr^(r, t) , 
An r (r, and An e (r, t) have been evaluated for 
r = 0.1/z, I/a, and 10/ll, and T h (r , r)=600 °C; for 



seven values of t, 3 ns ^ t ^ 3 ms; for ten values of 
r, r ^ r ^ 0.2 cm; and for ten values of time 3 ns 
^ t ^ 3s, where r, r, and t are all equally spaced in a 
logarithmic manner. The presentation in this paper 
of all these numerical results would be overwhelming. 
Table 10 contains a summary of these results for the 
case of platinum in Glass(B). The results for antimony 
and Al 2 3 (a) in Glass(B) and Glass(U) and for platinum 
in Glass(U), are similar to within a factor of ten. 



Table 6 

The maximum tensile stress as a function of r for fixed energy 
density E L = 20(J/cm 2 ) and pulse width t = 30 ns. The theoretical 
breaking stress of glass is about 6X 10 9 (N/m 2 ). The properties of 
Glass(B) are used in this table. The radius of the inclusion is r ; the 
maximum tensile stress is cr^n (max-tensile) and occurs at position r 
and at time t; the temperature at r and t is T h (r, t); the volume 
specific heat ratio is R; and the effective expansion coefficient is 
6a e ff. 



Inclusion: Pt(B), R = 2.97, 8a e 



-5.5X10- 6 (1/°C) 



logio(r ) 


logioW3) 


(r/r ) 


T h (r,t) 


(Too (max-tensile) 


(r in cm) 


(t in s) 




(°C) 


(N/m 2 )Xl0 9 


-2.0 


-8.0 


1.1 


0.0 


0.13 


-3.0 


-8.0 


1.1 


.0 


1.3 


-4.0 


-8.0 


1.3 


910. 


7.6 


-4.3 


-8.0 


1.4 


5200. 


9.0 


-5.0 


-8.0 


3.3 


370. 


2.6 


-6.0 


-9.0 


11. 


98. 


0.098 


-7.0 


-9.0 


110. 


0.80 


.00097 



Inclusion: Sb, # = 1.13, 6o; eff = -6.7 X 10- 6 (1/°C) 



logio(r ) 


log,oU/3) 


(r/r ) 


T h (r,4) 


doe (max-tensile) 


(r in cm) 


{t in s) 




(°C) 


(N/m 2 )Xl0 9 


-2.0 


-8.0 


1.1 


0.0 


0.27 


-3.0 


-8.0 


1.1 


.0 


2.7 


-4.0 


-8.0 


1.3 


2100. 


13. 


-4.3 


-8.0 


1.3 


4300. 


13. 


-5.0 


-9.0 


1.8 


1600. 


3.4 


-6.0 


-9.0 


11. 


100. 


0.097 


-7.0 


-9.0 


110. 


0.80 


.00097 



Inclusion: Al 2 3 (a), R = 2.74, da ef{ = -3A X 10- 6 (1/°C) 



logi (r„) 


log, (*/3) 


(r/r {) ) 


T h (r,t) 


(Too (max-tensile) 


(r in cm) 


(t in s) 




(°C) 


(N/m 2 )Xl0 9 


-2.0 


-8.0 


1.1 


0.0 


0.19 


-3.0 


-8.0 


1.1 


.0 


1.9 


-4.0 


-8.0 


1.2 


4600. 


11. 


-4.3 


-8.0 


1.4 


5600. 


13. 


-5.0 


-9.0 


1.5 


3500. 


3.6 


-6.0 


-9.0 


11. 


98. 


0.10 


-7.0 


-9.0 


110. 


0.80 


.00096 
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Table 7 



Table 8 



The maximum tensile stress as a function of r<> for fixed energy 
density E L = 20(J/cm 2 ) and pulse width r = 30 ns. The theoretical 
breaking stress of glass is about 6 X 10 9 (N/m 2 ). The properties for 
Glass(U) are used in this table. The radius of the inclusion is r : 
the maximum tensile stress is cr«« (max-tensile) and occurs at 
position r and at time t; the temperature at r and / is Tt, (r, t); the 
volume specific heat ratio is R; and the effective expansion coefficient 
is 8a t >r ( . 



Inclusion: Pt(U), # = 0.97, 8a e{f = - 5.1 X 10- 6 (1/°C) 



log,o(/o) 


log„,U/3) 


(r/r ) 


T h (r,t) 


cr^(max-tensile) 


(r in cm) 


it in s) 




(°C) 


(N/m 2 ) X 10 9 


-2.0 


-8.0 


1.1 


0.0 


0.18 


-3.0 


-8.0 


1.1 


.0 


1.8 


-4.0 


-8.0 


1.2 


2800. 


10. 


-4.3 


-8.0 


1.4 


2900. 


11. 


-5.0 


-9.0 


1.5 


1900. 


2.8 


-6.0 


-9.0 


8.6 


87. 


0.081 


-7.0 


-9.0 


80. 


1.3 


.00081 



Inclusion: Sb,fl = 0.49, 8a eff = - 6.7 X 10-«(1/°C) 



logio(ro) 


log,oW3) 


(r/r ) 


T„(r.t) 


o-^(max-tensile) 


(r () in cm) 


(/ in s) 




(°C) 


(N/m 2 ) X 10 9 


-2.0 


-8.0 


1.1 


0.0 


0.27 


-3.0 


-8.0 


1.1 


.0 


2.6 


-4.0 


-8.0 


1.2 


4600. 


11. 


-4.3 


-8.0 


1.2 


4300. 


11. 


-5.0 


-9.0 


1.7 


1200. 


2.7 


-6.0 


-9.0 


8.6 


89. 


0.079 


-7.0 


-9.0 


80. 


1.3 


.00081 



Inclusion: Al 2 0. { (a), R = 1.19, 8a eff = - 3.4 X I0- 6 (1/°C) 



logio(ro) 


log„,(//3) 


(r/r ) 


T,(r, /) 


o-M(max-tensile) 


(r in cm) 


U in s) 




(°C) 


(N/m 2 ) x 10 9 


-2.0 


-8.0 


1.1 


0.0 


0.19 


-3.0 


-8.0 


1.1 


.0 


1.9 


-4.0 


-8.0 


1.2 


2400. 


11. 


-4.3 


-8.0 


1.3 


6900. 


12. 


-5.0 


-9.0 


1.5 


1700. 


3.3 


-6.0 


-9.0 


8.6 


86.2 


0.083 


-7.0 


-9.0 


80. 


1.3 


.00081 



Maximum tensile stress of Pt spheres in a fictitious host as a func- 
tion of the thermal expansion coefficient of the host. All the prop- 
erties of the host are the same as those for Glass(B) except that the 
thermal expansion coefficient of the host an, varies from 6X10 -6 
(1/°C) to 16X10- 6 (1/°C). The effective expansion coefficient is 
8a eff ; the maximum tensile stress at a distance r and at time t is 
(Too (max-tensile): and the temperature of the host is F/,(r, t). The, 
laser pulse width is 7 = 30 ns and its energy density is Ei — 20(.)/cm 2 ). 
The radius of the sphere is r ( i = 5x 10~ 5 cm. The maximum stress 
occurs for times close to the pulse width, / — r. 



Oth 


8a e{{ 


r/r 


cree (max-tensile) 


T,Ar. t) 


(1/°C)X 10- fi 


(l/°C)x 10- 6 




(N/m 2 )X LO 9 


(X 10 ; '°C) 


6.0 


- 1.6 


1.4 


9.0 


5.2 


8.0 


-4.2 


1.4 


8.9 


5.2 


10. 


-6.7 


1.4 


8.8 


5.2 


12. 


-9.3 


1.5 


9.1 


2.1 


14. 


-12. 


1.5 


9.3 


2.1 


16. 


-14. 


1.5 


9.7 


2.1 



Table 9 

Maximum tensile stress of Pt spheres in a fictitious host as a func- 
tion of the thermal expansion coefficient of the host. All the prop- 
erties of the host are the same as those for Glass(U) except that the 
thermal expansion coefficient of the host an, varies from 4 X 10~ 6 
(1/°C) to 12X10- 6 (1/°C). The effective expansion coefficient is 
8a eff : the maximum tensile stress at a distance r and at time t is 
(Too (max-tensile); and the temperature of the host is T fl (r, t). The 
laser pulse width is r = 30 ns and its energy density is £/.= 20(J/cm 2 ). 
The radius of the sphere is r,, = 5X 10 -5 cm. The maximum tensile 
stress occurs for times close to the pulse width, t ~ r. 



Oth 


8a eft 


r/r 


(T (max-tensile) 


T„(r.t) 


(l/°C)x 10- 8 


(1/°C)X10- 6 




(N/m 2 )X 10" 


(X 10'* °C) 


4. 


1.2 


1.3 


12. 


7.9 


6. 


-1.3 


1.3 


11. 


7.9 


8. 


-3.9 


1.4 


10. 


2.9 


10. 


-6.4 


1.4 


11. 


2.9 


12. 


-9.0 


1.4 


11. 


2.9 



Table 10 contains the optical path length change 
AL and the stress birefringence parameters (Arc, 
and \nu) which occur when Pt is in Glass(B) for several 
values of pulse width t and inclusion size r . The 
results for Pt in Glass(U) are similar within a factor 
of about two to those predicted for Glass(B) and are 
not included here. Because the model assumes that the 
strain is proportional to stress, the energy densities 
El of the laser pulses are limited in table 10 to values 
such that the temperature of the glass never exceeds 
the strain point of the glass. That is, the energy 
densities are such that Th(ro, r) = 600 °C. In addition, 
the maximum tensile stress is always at least an order 
of mangitude less than the theoretical breaking point 
of the host 6 X 10 9 (N/m 2 ). Hence, excluding the 
complications which mierocracks might produce, one 
does not expect damage to occur when Tf,(r {) , r) 
^600°C. The quantities (AL), <An r ), and (An e ) 



denote representative values of the quantities AL(ri, 
r , t), Arc r (r, t), and Arc^(r, t) for those regions of 
r () ^ r ^ r 2 and t\ ^ t ^ to in which they obtain their 
largest values. The meaning of these symbols is 
illustrated below for one case from table 10. For ex- 
ample, the largest optical path length changes due to 
a platinum sphere of radius r = 1/x heated by a laser 
beam having a pulse width 3 X 10" 8 s and an energy 
density 0.33 (J/cm 2 ) are about 1.3 X 10 -8 cm and occur 
for times between 3/xs and 30/xs after the beginning 
of the laser pulse and for optical paths which approach 
within 1/x of the surface of the inclusion. Similarly, 
the largest values for the radial and tangentail changes 
of the index of refraction due to thermal birefringence 
are respectively 6.3 X 10~ 4 and 6.0X10 -4 and occur 
for the same radial distance and time as those for the 
optical path length changes. 
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Table 10a 

Optical path length change and stress birefringence for Pt inclusions in Glass(B). The radius of the inclusion is 
r (cm), the pulse width is t(s), the energy density is £Y(J/cm 2 ), and the optical path length change at time t=r and for a 
distance of closest approach r=ro is AL(r , r , r)(cm). The quantities (AL)(cm), (An r ), and (Arc^) denote respectively 
representative values for the largest optical path length change, the largest index change for radially polarized light, 
and the largest index change for tangentially polarized light which occur for a range of distances of closest approach 
r () to r 2 (cm) and for an interval of time t\ to £ 2 (s) after the beginning of the laser pulse. The tangential stress at r= r 
and t — T is — 2.16X 10 8 (N/m 2 ) and the surface temperature T h {r , t) at t = r is 600 °C in all cases. The radius is 
r«= 10~ 5 cm. 



logio(r/3) 


E L 


log,«(^/3) 


logio(*2/3) 


logio(r 2 ) 


AI(r () , r ,T) 


(M) 


<An,> 


(An„) 


(rin s) 


(J/cm 2 ) 


(t\ in s) 


(h'm s) 


(r 2 in cm) 


(cm) X 10- 6 


(cm) X 10- 6 


X10- 4 


X10- 4 


-3 


6035. 


-4 


-3 


-4 


0.24 


0.20 


2.1 


1.7 


-4 


605. 


-5 


-4 


-4 


.23 


.14 


2.0 


1.6 


-5 


61.1 


-6 


-5 


-4.7 


.19 


.14 


1.9 


1.5 


-6 


6.32 


-7 


-6 


-4.7 


.14 


.088 


5.8 


4.5 


-7 


0.706 


-7 


-6 


-4.7 


.094 


.033 


4.7 


3.2 


-8 


.113 


-7 


-6 


-4.7 


.047 


.0049 


0.022 


0.021 


-9 


.0446 


-7 


-6 


-4.7 


.025 


.0019 


.0080 


.0076 



Table 10b 

Optical path length change and stress birefringence for Pt inclusions in Glass(B). The notation of this table is 
the same as that for table 10a. The radius is r<> = 10 -4 cm. 



log,o(r/3) 


E L 


log, (*i/3) 


log, (* 2 /3) 


logio(r 2 ) 


AL(r»,r«,T) 


<AL) 


(An r ) 


<A^) 


(rin s) 


(J/cm 2 ) 


(t\ in s) 


(t-2 in s) 


(r 2 in cm) 


(cm) X 10- fi 


(cm) X 10-« 


x 10- 4 


x 10- 4 


-3 


611. 


-4 


-3 


-3.3 


1.9 


1.5 


8.1 


6.3 


-4 


63.2 


-5 


-4 


-3.3 


1.4 


1.0 


7.8 


6.0 


-5 


7.06 


-6 


-5 


-3.3 


0.95 


0.48 


6.8 


4.9 


-6 


1.13 


-6 


-5 


-3.7 


.47 


.049 


0.023 


0.22 


-7 


0.446 


-6 


-5 


-3.7 


.025 


.019 


.085 


.082 


-8 


.331 


-6 


-5 


-3.7 


.0 


.013 


.063 


.060 


-9 


.303 


-6 


-5 


-3.7 


.0 


.013 


.058 


.055 



Table 10c 

Optical path length change and stress birefringence for Pt inclusions in Glass(B). The notation of this table is 
the same as that for table 10a. The radius is r = 10 -3 cm. 



log,o(r/3) 


El 


log,oUi/3) 


log, o(£2/3) 


log, (r 2 ) 


AL(r„, r ,T) 


<AL) 


(An,) 


(kne) 


(rin s) 


(J/cm 2 ) 


{ti in s) 


(J 2 in s) 


(r 2 in cm) 


(cm) X 10-« 


(cm) X 10- 6 


X10- 4 


X10- 4 


-3 


70.6 


-4 


-3 


-2.52 


9.5 


6.3 


9.5 


7.1 


-4 


11.3 


-4 


-3 


-2.52 


4.7 


0.78 


6.9 


3.9 


-5 


4.46 


-4 


-3 


-2.52 


0.25 


.68 


1.8 


1.5 


-6 


3.31 


-4 


-3 


-2.52 


.0 


.52 


1.3 


1.1 


-7 


3.03 


-4 


-3 


-2.52 


.0 


.48 


1.2 


9.7 


-8 


2.95 


-4 


-3 


-2.52 


.0 


.47 


1.2 


9.4 


-9 


2.93 


-4 


-4 


-3.52 


.0 


.46 


1.2 


9.4 



The data from tables 6, 7, and 10 also may be used 
to estimate the lens effect due to heated regions of 
the host. Whenever the refractive index increases with 
temperature {dnnldTh) > 0, the heated region of the 
host surrounding the inclusion might focus the same 
laser pulse or a succeeding laser pulse. Computing 
the effective focal length for a spherical shell with a 
spatially varying refractive index is too complicated 
for an estimate of the effect. Instead the following 
crude approximation is considered. The focal length / 



of a spherical lens having a radius r and an uniform 
refractive index n\ is 

/= (ni/»o)r/[2{(7ii/iio)-lp], 

where n is the refractive index of the medium in 
which the spherical lens is imbedded. A measure 
of the expected importance of the spherical lens effect 
is the spatial extent of the deviation of the refractive 
index from its value before the laser pulse occurs, nn. 
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That is, the importance of the lens effect is measured 
in some manner by the quantity, 

1 [p 
(5ft) = lim - {/i/,(r, t) — n h }dr. 

P^ x P Jr 

A representative distance for the spatial extent is 

1 f p 
(r) = lim - r {/*,,(/% t) — n h }dr, 
p— p Jr 

where the above integrals are finite. The quantity 
[ft r (r, t) — ft/,] is proportional to the temperature 
7ft(r, 0; namely, [ft/,(r, - /i/,] ~ (dn h ldT h )T h (r, t) . 
When (5ft)— ^ft/,, then an estimate for the focal 
length is 

/~n«<r>/2<8*>» 

where 



and 



<r) = (dnuldTu) lim - f 'r^fr, *)<*r, 
p— pjr 

(r, r)dr. 



1 f" 
(8//) = (drih/dTh) lim - 7 1 /, 



Some typical numerical examples for platinum are as 
follows: when r = lp,, £/. = 20(J/cm 2 ), and r = 30 ns, 
then the maximum values for (r) and (5ft) occur at 
times t ~ 3jLts and are respectively (r) ~ 2jlc and 
(8ri) ~ 4 X lO- 3 . These values give an effective 
focal length of about / ~ 14 cm. But from tables 6 
and 7 the tensile stress at t = r = 30 ns exceeds the 
theoretical strength of the host before the lens effect 
becomes most important. Consider now another 
example for platinum in which (Tee (max-tensile) is 
less than the theoretical breaking stress of the glass 
host. When r = 1/x,, E L = 0.33 (J/cm 2 ), r = 30 ns, and 
T f ,(r {) , r) = 600 °C, then the maximum values of (r) 
and (8ft) occur at times about 3^s and are respectively 
(r) ~ 2/jl and (8ft) ~ 1.2 X 10" 4 . These values give 
an effective focal length of about /— 1.6 X 10 4 cm. 
These focal lengths are very long. 

Hence the model predicts that the lens effect arising 
from heated inclusions probably does not cause 
damage. Those cases for which the maximum tensile 
strength is less than the theoretical tensile stress of the 
glass have minimum effective focal lengths which are 
much greater than any dimension of neodymium doped 
glass elements used in present glass laser systems. In 
those cases for which the tensile stress exceeds the 
theoretical strength, the tensile stress probably causes 
damage before the lens effect could cause damage by 
heating another inclusion or by initiating an intrinsic 
damage mechanism such as self-focussing. 

Equations (34) and (35) and the results in tables 6 
and 7 show that Sb and the fictitious Al 2 0,3(a) are more 
likely to cause damage than Pt. However, this differ- 
ence is marginal. Additional computations for Pt 
spheres in fuzed silica, 8a eff ~ + 6.0 X 10~ 6 (1/°C); 
indicate that the maximum tensile stresses occur at 



r = r and t = r and that they are comparable to those 
stresses reported in tables 6 and 7. However, the range 
of r for which dee (max-tensile) exceeds the theoretical 
strength is much greater; namely, when 



io- 



r = 



10 6 cm, 



r = 30 ns, and E L = 20(Jlcm 2 ), then (Tee (max-tensile; 
Pt in silica) > 6 X 10 9 (N/m 2 ). Fused silica has a tensile 
strength comparable to glass. 

The maximum tensile stress as a function of the 
thermal conductivity Kn and the thermal expansion 
coefficient ah is studied also. These investigations are 
limited to Pt in Glass (B) and in Glass (U) for which the 
effective expansion coefficients are negative. It is found 
that increasing Kh in table 2 for Glass (B) from 0.008 
to 0.04 (W/cm °C) decreases the maximum tensile 
stress from 9.0 X 10 9 (N/m 2 ) to 4.0 X 10 9 (N/m 2 ). Simi- 
larly, increasing Kt, in table 2 for Glass (U) from 0.013 
to 0.04 (W/cm °C) decreases the maximum tensile 
stress from 10.6 X 10 9 (N/m 2 ) to 5.8 X 10 9 (N/m 2 ). 
Again, all the other properties of Glass (B) and of Glass 
(U) are kept the same except for the thermal expansion 
coefficient and the resulting maximum tensile stress 
as a function of ah for fixed r = 5 X 10~ 5 cm, r = 30 ns, 
and Ei = 20 J/cm 2 is reported in tables 8 and 9. Observe 
that the maximum tensile stress has a minimum at 
ah = ah(mm); namely, ^//(min, Glass (B)) ~ 10 X 10~ 6 
(1/°C) and a„(min, Glass (U)) ~ 8 X 10 6 (1/°C). Also, 
these tables indicate that the maximum tensile stress 
is a slowly varying function of the thermal expansion 
coefficient, a/,. Thus, for these cases, at, influences 
greatly the behavior near the interface because it 
appears in the effective expansion coefficient and the 
thermal conductivity influences greatly the behavior 
near the region of maximum tensile stress. 

Intuitive arguments exist to explain why the tempera- 
ture Tf, (r , i) in tables 6 and 7 has a maximum value 
at some r = r(max 7/,) for fixed E L and r and decreases 
for values of r greater than r(max T/, ) and less than 
r(msixTh). As (r /r( max Tf,)) becomes larger than one, 
the volume increases much faster than the surface 
area of the sphere. The inclusion receives less energy 
per unit volume. Therefore, the maximum surface 
temperature at the end of the pulse decreases. As 
(r /r(max Th)) becomes smaller than one, the equili- 
bration time (rg/ajO approaches zero. Then the surface 
temperature at the end of the pulse cannot deviate 
much from the equilibrium temperature due to the 
extremely short equilibration time. 

The above computations in table 10 suggest that 
examining incipient absorbing centers in laser glasses 
either by methods which employ the interference of 
two light rays (which experience different optical path 
length changes due to the local variation of the refrac- 
tive index near the inclusion) or by methods which 
employ thermal stress birefringence are promising. 
However, because the propagation of light in the pres- 
ence of a stress birefrigence which has spherical sym- 
metry is complex, the interpretation of birefringence 
data will be more tedious than data from an interfer- 
ence method. The feasibility of combining the inter- 
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ference method and holographic techniques has been 
demonstrated for 50 micron particles [2]. In addition 
the results of table 10 suggest that the use of laser 
pulses with pulse widths greater than a few micro- 
seconds may be more promising for the detection of 
small incipient absorbing centers than the use of nano- 
second laser pulses. The longer pulses produce spatial 
changes in the refractive index which extend over 
greater distances in the host and thereby increase the 
probability of detecting small inclusions. Another 
approach is to employ pulse widths and observation 
times which are less than the relaxation times of the 
host. Relaxation times for laser glasses between 800 
°C and 1100 °C are approximately between nano- 
seconds and microseconds. These are estimates of the 
time during which stress is proportional to strain in 
glasses which are not elastic for infinite time. Such 
short time observations will permit one to raise the 
inclusion surface temperature substantially above the 
present 600 °C reported here and still satisfy the 
assumptions of the present model. 
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Kahn for their reading of the manuscript. 



6. Appendix. Center-Host Interface 

A laser beam with an energy flux per unit area / 
(J/cm 2 s), with a pulse width r(s) and with a wavelength 
X/Xcm) impinges upon the inclusion. The inclusion 
exhibits an absorptance/* (A, T), a reflectance R(k, T), 
and a partial emissivity E(\, T), where X is the wave- 
length and T is the absolute temperature of the inclu- 
sion. The following relations among A, R, and E exist: 

(a) A-\-R = l; no energy is transmitted. 

(b) E(k, T)IA(K T)=e(K T); Kirchoffs law. 

The function e(X, T) is a universal function only of X 
and T and is independent of material and surface 
properties. 



(c) 



/; 

Jo 



e(X, T)d\ = SnT 4 ; Stefan-Boltzmann law. 



The constant S B is S B = 5.673 X 10~ 5 (erg/cm 2 s K 4 ). 

(d) X max r= 0.2897 cm K; Wien's displacement law. 
The wavelength for which e(X, T) has a maximum 
value for a fixed temperature T is X max . 

The values for the absorptance A(\, T) in table 4 
are valid for r () > X ~ 1.06 \x. The results of Mie scatter- 
ing theory [11] show that A(\, To) does not change 
appreciably for values of r () < X and r > y _1 where y 
is the absorption coefficient. 



The spherical inclusion has a complex index of 
refraction m c = n c — in' c and the Poynting vector for 
the incident radiation in the absorbing center is pro- 
portional to the factor exp (— y c (ro — r)). The absorp- 
tion coefficient y c is given by (47nic/A/J. Because a 
sphere with index m c is imbedded in the host with a 
real index mj, = rih ~ 1.52, the values for m c to use in 
the equations of ref. [11] are m c (med)=(/7i c //i/,). Simi- 
larly the wavelength in the host is X/, (med) = (\ijmi,) 
where \/,= 1.06 /t. The sphere intercepts in accordance 
with Mie scattering theory Q ex tirrll watts from the 
incident laser beam, independently of the polarization 
of the beam. There will be Qahs 77 ^ watts absorbed by 
the sphere and Qsca 77 ^ watts scattered in all direc- 
tions by the sphere. The conservation of energy gives 
@ext=(?abs+<2sca. The absorptance Q ahs = A(k, T; 
r () = small) replaces the absorptance A(k, T; r = large) 
quoted in table 4, whenever r becomes less than the 
wavelength X/(med). The variable X = (277r /X/Xmed)) 
is introduced and when X < 0.6 Mie scattering theory 
relates X and Q ext ; namely, 

Qext ^ Oabs « 2.68 X and Q sca < Q ahs . 

for platinum spheres in glass (rc/, = 1.52) and 

<2ext ~ Cabs ~ 0.992 X 

for antimony spheres in glass (tz/,= 1.52). When 
X=L06/z,Z = 9.0lX 10 4 cm- 1 r 
<2abs(P0 -24.1x104 cm- 1 r () 
<2abs(Sb) -8.93X10 4 cm- 1 r 



then 
and 



for values of r 5S 6.7 X 10~ 6 cm. Hence, A(k L , T) does 
not change by more than a factor of three for values of 
ro^yF 1 . 
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